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ABSTRACT 

«*>  ff 

detailed  analysis  of  the  nearly- frozen,  noneqiilllb- 

A 

rim  supersonic  expansion  of  a  dissociated  dlatcnlc  gas  around  a  sharp, 
tvo-dlnenslonal  comer,  assmlng  an  Inccnlng  stream  In  dissociation  equL- 
llhrlm.  G]he  analysis  Is  based  on  a  small  pertxirbatlon  technique  In  uhlch 
small  local  departures  from  frozen  flow  near  the  comer  are  treated.  Sie 
anal^ls  Is  necessarily  restricted  In  lateral  extent  but  covers  the  entire 
angular  range  of  the  flow  field  and,  therefore,  brings  ovtt  the  nonlinear 
aspects  of  the  nonequillbrlm  behavior.  An  exact  analysis  Is  presented 
(in  similitude  form  for  a  diatomic  gas)  In  which  numerical  procedures  are 
required  to  effect  the  solution.  However,  in  the  exact  solution,  the 
expression  for  the  atomic  species  mass  fraction  has  been  placed  in  closed 
form  euid.  In  conjxmctlon  with  appropriate  approximations.  Is  used  to  de¬ 
termine  the  theimochemlcal  field  In  an  additional  solution  (applicable  to 
hypersonic  flow)  which  is  approximate  but  entirely  in  closed  fom..  borne 
nunerical  results  for  the  first-order  nonequilibriun  flow  field  are  pre- 
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NCMEWCLATURE 


Sjnnbol 

a. 

A 

6 

<2 

r 
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Speed  of  Bovmd 
Constant  [see  Eq*  (2-13)] 

Constant  [see  Bq:.  (2-13)] 

Constant  pressure  specific  heat 
Frozen  constant  pressure  specific  heat 
Chemical  potential  (Glhbs  function) 
Function  defined  in  Eq.  (2-13) 


if 

h 

i^m 

K 

Kjr 

Kh 

M 


Specific  enthalpy 

Specific  dissociation  energy 

Integral  expression  defined  in  Eq,.  (3-3^) 

Integral  expression  defined  in  Ei^.  (3-3?) 

Recombination  rate  constant 

I^raneter  defined  in  E(|.  (2-2l) 

Mach  nvmber 
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NOMENCLATURE 


Symbol 

f 

P 

r 

R 

P. 

5 


VresBwre 

fife 

Radial  coordinate 
Gas  constant 
Uilversal  gas  constant 
Specific  entroiy 


J  5/R- 


T“  Temperature 

7]^  Characteristic  dissociation  temperature 

(j  Radial  velocity 

ff  y/vi,. 

V  Tangential  velocity 


V 


oc 


y 

y 

A 

y 

r' 


v/v,. 

Incoming  total  velocity 
Atomic  mass  fraction 
Actual  specific  heat  ratio 
frozen  specific  heat  ratio 
Effective  specific  heat  ratio 
Characteristic  length  [see  E^.  (2-22)] 


NCKESCUTURE 


Symbol 

•f  Vortlcity 

f 


[(5.>/yc7«-/)] 


*k 


■d- 

f 


Angular  coordinate 
Mach  angle 
FLOW  turning  angle 
Density 


O’ 

T 

r 

6J 


e/f- 

T/T^ 

Reduced  angular  coordinate 
Beccniblnatlai  rate  temperature  exponent 
Net  rate  of  atomic  species  mass  production 


Sxxbscrlpts 

A  Refers  to  atcmic  species 

f  Refers  to  frozen  or  zerotb  order  field 

X  Index  (generally  refers  to  the  chemical  species) 

M  Refers  to  molecular  species 

/  Refers  to  first  order 


Subscrl'pts 


I 

Denotes 

con(3itlons 

H 

Denotes 

conditions 

eo 

Refers  to  Incoming 

at  Initial  frozen  cliaracteristic 
at  final  characteristic 
stream 


1.  lUTRODUCTION 


problem  of  the  centered  expansion  of  ste^idy^  stqpersonlc,  reacting  gas 
flow  around  a  sharp  comer  (Prandtl -Meyer  flow)  is  of  practical  interest 
as  a  simple  case  of  nonequilibrim  flows  that  are  encountered  on  bodies 
in  high  altitiiie,  hypersonic  flig^ht  (Refs.  1,2).  Moreover,  the  stvuSy  of 
this  problem  constitvttes  a  fundamental  research  tool  with  which  to  inves¬ 
tigate  the  effects  of  homogeneous  reactions  on  supersonic  conical  flow 
fields.  Complete  solutions  for  a  dissociated  gas  mixture  have  been  given 
in  the  case  of  an  equilibriun  expansion,  \diere  the  dissociation  and  recom¬ 
bination  rates  are  equal  and  much  greater  than  the  convection  rate  throvi^- 
ovrt  the  expansion  fan  (Refs.  3,  ^,5).  presence  of  a  chemical  reaction 
in  local  eqviilibrium  does  not  change  the  basic  self-similar  nature  of  the 
solution  that  is  observed  in  a  perfect  gas,  since  the  flow  remains  isen- 
tropic  and  irrotational  throu^out  with  the  characteristics  emanating  as 
strai^t  lines  from  the  comer.  However,  when  the  incoming  flow  velocity 
is  sufficiently  hl^  and  the  density  sufficiently  low,  the  expsmsion  cam 
force  the  local  reaction  rates  out  of  equilibrium.  The  finite  dlssoclatlon- 
reccmbination  rates  introduce  a  nonlsentropic,  rotational  effect  into  the 
problem,  as  well  as  a  length  scale  proportional  to  radius  that  destroys 
the  conicity  of  the  flow.  With  increasing  radius,  there  is  a  recombination- 
dominated  (nonequilibrium)  transition  from  completely  frozen  flow  near  the 
comer  to  am  asymptotically-approached  region  of  equilibrium  flow  far  from 
the  corner.  When  the  incoming  velocity  is  very  hl^  (auid/or  density  very 
low),  the  entire  expansion  may  become  chemically  frozen,  in  which  case  the 
flow  aigain  reverts  to  an  isentropic,  irrotatlonaO.,  self-simllaur  behavior. 

Ihe  theoretical  description  of  a  nonequillbrlua  centered  expauislon  Involves 
a  set  of  coupled,  nonlinear,  partial  differential  equations,  therefore, 
in  general,  solutloxis  must  be  obtained  by  a  nunerical  technique  based  on 


the  method  of  characteristics  for  reacting  gas  flows  (Refs.  3,  7>  8). 
However,  valuable  insist  into  the  nonequi librium  effects  can  be  gained 
by  considering  small  departures  from  a  self-similar  solution  by  either 
(a)  treating  weak  expansions,  with  small  departures  frcm  conditions  In 
the  incoming  gas  stream,  or  (b)  treating  small  local  departures  fran  the 
frozen  flow  near  the  comer.  Method  (a)  features  greater  simplicity 
because  the  flow  may  be  regarded  as  irrotational  and  the  restilting  linear¬ 
ized  equations  possess  constant  coefficients.  Ihis  case  has  been  treated 
by  several  Investigators  for  both  expansion  and  canpression  flows  (Refs. 

9,  lO).  Ihe  method  is  restricted  to  small  turning  angles  but  covers  the 
entire  lateral  extent  of  the  field.  Method  (b)  is  mathematically  more 
ccmplex  since  the  coefficients  in  the  linearized  equations  governing  the 
nonequilibritm  perturbations  are  complicated  functions  of  the  zeroth  order 
self -similar  solvrtion.  Moreover,  the  perturbations  are  not  irrotational. 
Nevertheless,  this  approach  is  unrestricted  as  to  flow  turning  angle  and 
brings  out  the  nonlinear  aspects  of  the  nonequllibrlun  bdiavlor.  This 
report  presents  a  detailed  analysis  of  a  nearly- frozen,  nonequllibrlun 
supersonic  expansion  of  a  dissociated  diatomic  gas  around  a  sharp,  two- 
dimensional  comer  for  the  case  of  an  equllibrlw  dissociated  incoming 
stream.  Ihe  present  Investigation  cor^stitutes  a  further  extension  of 
Kapolitano's  analysis  (Ref.  8)  in  that  the  entire  angular  range  of  the 
flow  field  is  treated,  sijecific  numerical  results  are  given,  and  seme 
approximate  and  exact  closed  form  solutions  are  developed. 

2.  GOVERNING  RELATIONS 


2.1  Basic  Equations 

Consider  the  steady,  adiabatic,  supersonic  expansion  of  a  dissociated 
diatomic  gas  around  a  sharp  two-dimensional  comer,  as  shown  in  Tig,  1. 
Viscosity,  heat  conduction  and  diffusion  effects  sure  neglected,  smd  the 
flow  properties  Immediately  ahead  of  the  first  expansion  fan  chsuracteristlc 
sure  stssxaned  to  be  xml  form,  nie  governing  equations  for  such  a  reacting 
gas  flow  ax’e  well-known  (see,  for  example.  Ref.  u).  In  the  polar  coor¬ 
dinate  system  of  Pig.  1,  they  are  expressed  as  follows. 
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Continuity 


.^-.0 


<2-l) 


Aton  Mass  Conservation 


Radial  Monentm 


Tangential  Mcmentm 

/s/i.dv  .  V  dv  .  uv  ^  _  /  a-p 

+ -’T-SS  ♦ -Y^;  -  -  r 

Energy 

/?  'h  s  CONST 


(2-2) 


(2-3) 


(2-4) 


(2-5) 


Theimal  Eqtuation  of  State 


-P  ^  ef^^d'hocyr  (2-6) 

Ihe  specific  enthalpy  h  is  given  as  a  function  of  oC  and  T  by  the 
caloric  equation  of  state: 


(2-T) 
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vhere  C  ^  denotes  the  atonic  or  molecular  specific  heat  and  is 
the  dissociation  energy  of  the  gas.  In  the  subsequent  analysis,  the  rota¬ 
tional  degrees  of  freedom  of  the  molecule  are  assumed  to  he  fully  excited, 
and  the  effects  of  electronic  excitation  and  vibrational  nonequilibrium 
on  the  internal  energy  are  neglected  in  comparison  to  the  dissociation 
energy  oi  hp.  ^lus  Cp^  “  ®  function  of  T  ranging  be¬ 
tween  (no  vibrational  excitation)  and  -fr  (completely  excited 

vibration).  To  further  simplify  the  analysis,  Cp^  shall  be  regarded  as 
a  constant  (l.e.,  vibration  is  "frozen"  throuf^out  the  expansion).  Hence 
Bq.  (2-7)  becomes 

/)  =  C”p  T  +  ,  (2-8) 

where  Cp  =  is  the  frozen  specific  heat  and 

i 

Die  entropy  chaafje  along  a  streamline  for  reacting  diatomic  gas  flow  is 
(Refs.  6,  U): 

’  (2-9) 

■vdiere  ^  -  T  5  ^  ^  denotes  the  "chemical  potential"  (or 

M  A,  ^ 

Gibbs  function)  for  the  atom  or  molecule  (expressions  for 
and  5  in  tenns  of  T,  -p  and  OC  are  given  in  Appendix  A).  The  follow¬ 
ing  equivalent  form  of  Eq.  (2-9),  namely  Crocco's  Iheorem  for  a  reacting 
gas,  may  also  be  obtained  from  (2-9)  and  the  momentum  equations: 

(2.10A) 

r-H-uf  ,  (2-ioB) 

I 
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\Aiere  -f-  vorticlty  in  cylindrlcaa  coordinates.  Both 

Eqs.  (2-9)  and  (2-10)  are  very  useful  in  obtaining  a  physical  appreciation 
of  the  nonequillbrivsn  dissociation  effects  on  an  Invlscid  flow. 

2.2  Reaction  Rate 

Ihe  dlssoclatlon-recanblnatlon  reaction  for  a  diatonic  gas  is 


x  +  M  ^  24  4.  X  ,  (2-11) 

\diere  M,  A  denote  a  molecular  or  atomic  particle  and  the  third  body  X 
may  be  either  of  these  two.  For  such  a  reaction,  the  net  volwietrlc  rate 
of  atomic  species  mass  production  can  be  written  as 

s  -  ,  (2-02) 

where 

»  /  -r-  “Vr 

=  7?^- r/- oc)(^j  ^  e  (2-13) 

Here,  and  cJ  represent  the  recombination  rate  constant  and  tempera¬ 
ture  exponent,  respectively,  A  and  B  are  constants  (Appendix  A),  and 
is  a  characteristic  dissociation  tenqperature .  Ihe  function 
G’fo4,^,7’^  has  the  Important  property  of  vanishing  identically  when 
themochemlcal  equllibrlun  prevails.  !Ihus  when  reaction  (2-ll)  is  in 
local  equilibrium  in  the  flow,  the  equation  Gr  ^  O  defines  the  correspond¬ 
ing  solution  to  (2-2)  for^  . 

2.3  Speed  of  Sound 

Ihe  speed  at  which  infinitesimally  small  disturbances  propagate  through  a 


5 


dissociated  diatomic  gas  can  be  written  (Ref.  6)  as 


« 

where  J  “  0  when  the  gas  is  not  in  chemical  equilibriun  and  J  -  /  when 
the  gas  remains  in  equilibriun  ipdien  subjected  to  an  infinitesimal  distur¬ 
bance.  Here,  y  is  the  actual  specific  heat  ratio  of  the  dissociated  gas, 
defined  by 


y  s _ _  , 


and  y  is  the  frozen  specific  heat  ratio 


(2-15) 


Ys  (y)  a  — _  a  ^ . 

-r.caNST  I^(^-Km')<X 


(2-16) 


According  to  (2-l4),  the  velocity  normal  to  a  Mach  wave  equals  the  local 
"frozen”  speed  of  sound  (based  on  y )  in  any  gas  with  a  finite  reaction 
rate.  However,  in  a  ccmpletely  eq\iilibriuii-diBsoclated  gas,  the  speed  of 
sound  based  on  y  (y<  y  £.  y)  must  be  used. 


2.U  Boundary  Conditions 

If  the  incoming  gas  is  uniformly  dissociated,  so  that  the  Initial  expansion 
characteristic  is  strale^t,  it  must  be  either  chemically  frozen  (after 
having  previously  undergone  a  hic^ly  nonequilibriun  flow  expanslcai)  or  in 
thermochemical  equlllbrlvm.  In  the  former  case,  oc^  is  an  independent 
variable  with  while  in  the  latter  case,  is  & 

function  of  and  7^  given  by  s  O  •  Since  any  expan¬ 

sion  of  a  chemically  ffozen  stream  must  remain  frozen,  we  shall  consider 
only  an  equilibriun  fJree  stream.  Consequently,  th«>  following  conditlcxis 


define  our  problon  as  a  supersonic  expansion  around  a  Sharp  comer  (fig. 
l).  At  the  initial  expansion  fan  characteristic  I, 

■e-r'A 


"f*  ■  />oo  j  p  «Pa»  y  T  *  7»  t  ^  *  OC  ca  1 

U  =  Cos  ,  V  "^i  ) 

At  ^  s  - 1/",  the  tangential  velocity  ccenponent  must  vanish : 


(2-18) 


(2-19) 


The  flow  along  the  wall  for  a  nonequilihrlum  expansion  will  he  a  fimction 
of  r*  since  a  relaxation  toward  equilibrim  occurs  with  increasing  stream¬ 
line  distance. 


2.5  Ropd^^nensional  Equations 

It  proves  convenient  to  deal  with  the  foregoing  equations  in  an  appropri¬ 
ately  nondlmenslonallzed  form.  Therefore,  the  following  new  Independent 
variables  are  Introduced: 


fi  =  i>lfm  ,  cr=  e/e« ,  T -  7-/7i 

y-U/Viv,  '  >Jr- 


Furthermore,  we  define  the  parameters 


r..r./r.,K,  =  -^f5^,  K,=  W^. 


(2-20) 


) 


(2-21) 


and  a  characteristic  (flow  length/reccoiblnation  length)  parameter 


-r  ^rr 


(2-22) 


Ohen  Bq.8.  (2-l)  -  (2-6),  using  (2-8),  (2-32)  and  (2-13),  are  transformed 
into  the  following  set  of  six  equations  for  the  six  unknowns  od ,  P  , o*  , 
T,  U  and  V  ; 


^-0 


(2-23) 


(2-25) 


ap 


(2-26) 


U V V* /)  7^[rp  +  ('f-KH)Tj  (2-27) 


Ihe  attending  entropy  relations  become 


(2-28) 


(2-29) 


or,  In  terns  of  Croceo's  ISieoren, 


’'■Ir  * 

vihere 


(2-30A) 


(2-30B) 


Note  that,  according  to  Appendix  A,  ^  vritten  in  tema  of 

the  reaction  function  G-  as 


rSrf K— •  ^*i]  >  (2-31) 

vhlch  properly  vanishes  at  equilibrlvnn  {G~0)»  Finally,  the  boundary  con¬ 
ditions  transform  to 

^  6^) : 


Pscr=  T  =  ( 


oc  «  oc. 


/-o<j  =  e 


-■«h> 


u  * 


(2-32) 
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■P-  =  -  2^  : 


(2-33) 


The  foregoing  equations  are  in  a  universal  nondlnenslonal  fozm  which  is 
very  convenient  for  parametric  study  of  nonequilibriun  expansions  regard¬ 
less  of  the  method  of  solution.  Furthermore,  it  is  evident  that  P  ,“C , 

0^ ,  CT  ,  U  and  V  at  a  given  and  -xr  are  the  same  at  each  point  -O', 
r'T  for  all  combinations  of  inccmlng  gas  if  the  following  parameters  are 
held  constant:  Kp.  That  is,  the  following  similitude  law 

governs  the  flow: 

oc, o-,U,V‘  f  )  (2-34) 

This  simllltiide  is  in  many  respects  analogous  to  that  demonstrated  in  Ref. 
22  for  nonequilibrium,  small  disturbance  hypersonic  flows. 

It  is  seen  that  the  flow  described  by  the  foaregolng  equations  cannot  be 
conical  in  the  presence  of  a  finite  dissoclation-z^ccmblnation  rate,  since 
the  parameter  P  is  a  function  of  T*  .  Purthennore,  "vdien  H and 
G^O,  Bqs.  (2-24),  (2-29),  and  (2-30)  show  that  the  degree  of  dissocia¬ 
tion,  entropy  and  vorticlty  vary  along  streamlines  as  the  flow  expands 
around  the  comer.  The  effects  of  nonequllibritm  chemical  reaction  Involve 
a  transition  between  frozen  flow  at  the  comer  iP»0  O  ,oc  -oC«)  and 
equlllbrlvm  flow  far  from  the  comer  (r-^co,^  =  O  ,ac«0cf^)  — )  across 
which  the  entropy  Increases  and  oc  decreases  with  7*  .  The  radial  expeuise 
of  this  transition  region  is  governed  by  the  magnitude  oi  P*  .  Vhen 
r  «  /  ,  the  radial  sceae  vahishes  frcm  the  problem  and  the  entire  expan¬ 
sion  beccmes  a  chemicaTJy  frozen,  conical  flow.  As  P*  is  Increased,  the 
otrter  edge  of  the  nonequillbriiin  z<»e  moves  inward.  When  T  »  /  ,  the 
entire  expansion  (with  the  exception  of  a  very  small  region  near  the 
comer)  is  in  equillbrlun  and  hence  is  again  self -similar .  In  this  case, 
of  course,  oc  decrecuses  throu^  the  expeuislon  according  to  the  rela¬ 
tions  G  "  O  and 
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(2-35) 


Ass.  - 

ct-d- 


oc 


JjL 

■- 


Die  presence  of  chemical  nonequilibrium  also  induces  a  curvature  of  the 
expansion  fan  characteristics.  Ilhls  cvirvature  may  be  observed  by  noting 
that  departures  from  frozen  flow  must  Intrclace  a  reaction  rate  effect 
which  is  recombination-dominated  and  which  therefore  tends  to  reduce  , 
luid  increase  TT  and  P  !Ihus  for  a  given  locsQ.  turning  angle ,  the  Mach 
number  Increases  with  radius  and  the  expansion  fan  characteristics  curve 
iqward  and  away  from  the  frozen  flow  characteristics.  At  a  sufficiently 
large  distance  from  the  comer,  the  characteristics  become  asjanptotlcally 
parallel  to  (but  not  coincident  with)  the  characteristics  for  a  completely 
equilibrivm  expansion  (Fig.  l).  This  lack  of  coincidence,  which  is  a  mani- 

— i 

festation  of  the  radial  increase  in  entropy,  vanishes  when  /  . 

3.  MEARLY-FROZEK  NONEQUIUBRIH!  EXPANSIONS 
3.1  Perturbation  Method 


We  now  consider  small  departures  from  a  completely  frozen  expansion  due 
to  nonequilibrlvm  dissociation  effects.  Since  a  nearly  frozen  flow  is 
characterized  by  small  values  of  P  *  P  T*  ,  we  expand  each  of  the  depen¬ 
dent  variables  in  a  series  of  ascending  powers  of  P  : 

ocC^,t)  =oc^C-^)^-pa,C^)  .#.p*oCi  -  - 

*rp  c^)  ^  • 


Vtien  a  dissociated  gas  is  expanded,  the  nonequlUbrlun  effects  on  resultant 
velocity  and  flow  density  are  less  significant  than  the  corresponding 
effects  on  temperature  and  pressure. 
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etc.,  where  subscript  F  denotes  the  local  frozen  expansion  solution, 
subscript  f  denotes  the  local  first  order  perturbations  due  to  nonequi¬ 
librium  effects,  etc.  Substituting  into  Eqs.  (2-23)  -  (2-28),  equating 
coefficients  of  like  powers  of  P  equal  to  zero  and  taking  K^*/  ,  we 
obtain  a  set  of  nonlinear  ordinary  differential  equations  for  the  zeroth 
order  (frozen)  flow  field  and  a  set  of  linear,  nonhcmogeneotus,  ordinary 
differential  equations  with  variable  coefficients  for  each  of  the  first, 
second,  etc.,  order  perturbations.  Confining  our  attention  to  first-order 
nonequllibrixxn  effects,  we  have  the  following  two  sets  of  equations. 

Zeroth  Order  (Trozen)  Flow 


(3-1) 

oCpC*©-)  S  CONST  :  C<a> 

(3-2) 

(3-3) 

CONSrr»L  ) 

(3-4) 

0  -  Tp)  =  (  Up  +  V^-0  (3-5) 

Pp  =  o-p  Tp.  (3-6) 
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First -Order  Flow 


Z(a;  £7p  +  ov  ^  (cTp  v<  +  cr,  )  =  0 


(3-T) 


=  -2/’/Tr^/5^  e^'J 

sZ _ _-L _ iy 


(3-8) 


r 


(3-9) 


(3-10) 

(^)T,y-/Si±^|iS^7«,^5,Mi(?fir*{?i7;=0  (3-u) 


i  =  -<x,  ^ 

Pp  op  TT^  /+0^«»  ”T*p 


Ihe  corresponding  entropy  and  vortlclty  perturbation  equations  are 


(3-12) 


U^J^-hV  = 
^F4f,  •»■  «V 


(3-13) 


p2  )Kpt/  e  7f  fft  7 


and 


(3-l^A) 


(3-15) 

Ihe  boundary  conditions  on  the  frozen  and  first-order  solutions  are 


m  p  ^  -  f—  —  1  ri  CMi“0.  ^  I 

.fr- ‘v-Of-',  V,-  >  »P 

(3-16) 

/;  =  T,  =<r,=ec,=  q  =  !^  .0 

(3-17) 

•^  =  -  ; 

II 

-<l 

H 

O 

(3-l8) 
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3.2  Frozen  Ebcpanslon  Solution 

Eqs.  (3-l)  -  (3-6)  describe  the  self -similar  expansion  of  a  theimally  and 
calorically-perfect  gas  with  a  specific  heat  ratio  3hus  the 

effects  of  frozen  dissociation  on  the  expansion  are  bounded  by  the  solu¬ 
tions  for  an  undissociated  perfect  gas  (^s  z  1.4)  and  a  pure  non- 
atcmic  gas  ( 5^#  *  ^  =  1.66?)  for  a  given  and  Since  the  details 
of  the  frozen  flow  field  are  needed  to  evaluate  the  coefficients  and  non- 
hcmogeneous  terns  in  the  first-order  perturbation  equations,  a  brief 
review  of  the  zeroth-order  solution  will  now  be  given. 

Ihe  analysis  is  naturally  divided  into  three  distinct  regions:  A,  B,  and 
C  (FLg.  l).  In  region  A  ahead  of  the  first  characteristic,  the  resultant 
velocity  and  theznodynamlc  properties  are  uniform  \Aille  the  tangential 
and  radial  velocity  c<xq>onents  are  given  by  pure  geometry.  Ihus,  when 


(3-19) 

Up  s  Cos  >  Vp  -  Sir4-0~ 

In  region  B  within  the  expansion  fan,  the  flow  is  determined  by  Eqs.  (3-l) 
throu^  (3-6)  and  (3-19)*  Combining  (3-l)#  (3-3)  and  (3-6),  we  obtain 


(3-20) 


which  Indicates  that  the  tangential  velocity  conponent  Is  equal  to  the 
local  speed  of  sound  In  region  B.  Substituting  the  Irrotatlonallty  con¬ 
dition  (3-3)  and  (3-20)  Into  the  energy  equation  (3-^)>  the  following 
first-order  differential  equation  for  Is  obtained: 


(3-fil) 


\diere 


Taking  Into  account  the  initial  houndary  condition  on  Upt  Eq.  (3-21) 
Integrates  to 


(3-22) 


vhere  the  choice  of  sign  in  is  determined  as  shown  helow.  Ihus,  using 
Bq..  (3-3)  and  the  "boundary  condition  on  Vp  ,  we  obtain 

Denoting  %  at 


Bqs.  (3-22)  and  (3-23)  can  be  written  as 

.  t  J 

n  -  (^m ^.0  ^luy 
Sin  Tr 


(3-24) 


J-CaiJr 

Mm  Cos  r. 


(3-25) 


According  to  (3**20);  the  frozen  temperature  dlstrlhutlon  In  the  expansion 
fan  becomes 


(3-26) 


and,  therefore,  since  the  flew  is  Isentropic, 


-  ( Cosy 
\  Cos 


(3-2T) 


(3-88) 


^e  Mach  nunber  distribution  in  the  fan  now  may  be  written  as 


+  /  (3-29) 

Ihe  solution  is  completed  by  observing  that  the  sign  on  -  -©i.  in  the 
foregoing  definition  of  mxist  be  such  that  ,  Pp  ,  <  /  (and 

M  >M^)  for-©-  <*6^*  Since forM^<oo,  BVis.  (3-26)  -  (3-29)  re¬ 
quire  and  hence  a  negative  sign  on-»~-6y.  i.e.. 

z  "i  I 


The  solution  in  the  third  region,  C,  the  unifom  flow  field  downstream  of 
the  last  frozen  expansion  fan  characteristic  which  is  parallel  to  the  wall, 
again  becomes  quite  simple.  The  theznodynamic  properties  and  resultant 


flow  velocity  are  equal  to  their  values  at  and  the  radial  and 

tangential  velocity  components  are  given  by 


Si«(e-+vv) 


(3-31) 


Using  (3-29)3  calculated  as  a  function  of  Mo*  and  ^  by  the 

relation 


Tau  Yy*>^[ 7 ,  (3.32) 

where  i®  "th®  Prandtl-Meyer  turning  angle  function  (Ref.  13)  and 

where  the  choice  of  sign  rests  on  obtaining  the  proper  iqward-running 
characteristic  (discarding  the  downward-rtinnlng  member  In  this  problem). 

3.3  First-Order  Solution 

Bq.  (3-8)  is  decoupled  from  the  remaining  first-order  equations;  a  formal 
solution  for  <K ,  can  therefore  be  easily  obtedned.  Switching  to  the 
independent  variable  ,  inserting  the  proper  frozen  flow  functions  for 
the  coefficients,  and  using  the  appropriate  integrating  factor  (3-8)  yields 
the  following  solution  in  region  B: 


CAiYjCosir)  IP,  r,  '  (CosY)  &r  jy 


(3-33) 


With  OC,(y^'^xO  •  Using  (2-13),  (3-26)  and  (3-27)3  and  noting  that 


18 


$ 


K^0-o(,) 


(3-33)  can  be  written  as 


f 


.w 


(3-34) 


where 


i 


(ccB^o  ar 


(3-35A) 


^W^jiSicy)  cosV,, Sec*y)<ky 


(3-35B) 


represent  contributions  Arcoi  the  reccniblnation  and  dissociation  rates, 
respectively,  the  fonoer  being  the  predominant  meehaxiism  since  oc  de¬ 
creases  with  T*('o(,£o)  .  Bq.  (3-35)  correctly  predicts  that  Oc,  vanishes 
when  the  incoming  stream  is  xmdissociated  (oi^»o). 

nie  foregoing  quadratures  mvmt  be  carried  out  ninerleally  for  arbitrary 
values  of  y  ,  and  A  ,  However,  when  ^^Zco  -  5  is  an  Integer, 
(3-35A)  may  be  evaluated  in  closed  fozm  directly  from  available  Integral 
tables  (Ref.  l4).  Moreover,  for  4  equal  to  an  even,  positive 

integer,  a  closed  fonn  representation  of  (3-35B)  in  tenns  of  the  error 
function  can  be  obtained  (details  are  given  in  Appendix  B)  .*  Fortunately, 


# 

Hie  authors  are  Indebted  to  Mr.  T.  Okanura  of  the  Research  Ghro\qp,  Missile 
Aero/TheimiodTnaBilcs  Section, for  pointing  out  this  representation. 
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the  conditions  for  closed  fom  integrabllity  lie  within  the  realjn  of  prac¬ 
tical  interest,  since  cj  is  usually  negative  {tosi  -1.5)  for  diatomic 
gases.  Therefore,  a  nmerical  evaluation  of  (3-35)  for  arbitrary  values 
of  ">1  ,  U)  and  A  can  be  avoided  by  a  parametric  study  of  the  closed  fom 
soltrbions  and  subsequent  interpolation. 


In  region  C  downstream  of  the  last  frozen  expansion  fan  characteristic 
the  solution  to  (3-8)  is  8trai^tforwca*d.  Here,  the  nonhomo- 
geneous  reaction  rate  tern  is  constant  and  U,- ,  are  given  by  (3-3l). 
The  resulting  integration,  subject  to  the  initial  condition 
can  be  carried  out  in  closed  fom  and  yields 


oc,s(oc.)  - 

^  S|m(^4>VV) 


(3-36) 


'Ij 


Sin(^  *yOl^oT 


This  solution  describes  a  reccmblnatlon-danlnated  relaxation  downstream  of 
the  final  frozen  expansion  fan  characteristic  in  which  the  Influence  of  the 
previous  nonequillbriun  flow  history  within  region  B  continually  diminishes 
as  decreases.  At  the  wall,  oC|  becomes  completely  Independent  of 

the  flow  in  region  B  and  nay  be  calculated  directly  from  the  frozen  solution. 

The  remaining  first  order  perturbations  nay  now  be  obtained  by  solving  two 
additional  differential  eqxiatlons.  Napoli tano,  in  his  discussion  of  the 
problem,  chose  to  evalmte  the  first-order  entropy  and  radial  velocity  by 
qxiadratures  and  then  determine  the  other  flow  field  variables  algebraically. 
This  approach  is  suggested  naturally  when  one  observes  that  (3-13)  Is 
uncoxqsled  from  the  remaining  eqwtlons  by  the  foregoing  solution  for  of , . 

In  practice,  however,  the  advantage  of  this  procedure  is  somewhat  diminished 

I 
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since  the  resulting  integration  must  still  be  dcme  nunerlcally.  Moreover^ 
as  shall  be  shown,  it  is  more  firultful  to  proceed  Instead  in  tezms  of 
the  vortlcity  field,  particularly  for  the  purpose  of  developing  approxi¬ 
mate,  closed  fonn  solutions,  ^erefore,  Eiis.  (3-13)  and  (3-1^)  will 
henceforth  be  discarded. 


Convenient  relations  governing  the  first-order  vortlcity  perturbation 


may  be  easily  derived.  First,  ve  eliminate  between  (3-6)  and 

(3-10),  using  (3-7)j  to  yield 


(3-37) 


A  second  vortlcity  relation  is  also  obtained  by  ccnblnlng  (3-9);  (3-ll) 
and  (3-32): 


Qiese  two  relations  are  equivalent  to  the  two  mcmentm  equations  (3-9)  and 
(3-10).  Row  within  the  frozen  ftm  region  B,  where  (3-37)  and 

(3-36)  may  be  combined  so  as  to  eliminate  the  coninon  factor  containing 
o;  ,  Uf  and  V|  ;  then  using  (3-l)  and  (3-3);  one  obtains  a  linear  first- 
order  differential  equation  for  that  readily  Integrates  to  the  solu¬ 
tion: 


2  ,  (7s>di  y 


(3-39) 


with  Of,  fv)  known  and  the  borjndary  condition  In  region  C 

downstreain  of  the  last  frozen  fan  characteristic,  the  solution  for  ^ 
becomes  quite  sluqple,  since  by  Eq.  (3-37) 


=  CONST  = 

lAen  ^^*0  •  Ihtu  the  nonequlllbrliiii  effect  within  region  B  xeeults  in 
a  flow  which  emerges  with  a  fixed  rate  of  rotation  proportional  to  T 
throu^out  region  C. 

With  ^  known,  the  differential  equation  governing  the  radial  wloclty 
pertiu-bation  Jj,  may  now  be  determined.  Die  expression  f,  = 
and  Bq.  (3-38)  are  differentiated  to  obtain  relations  for  and  , 

respectively.  These  relations  along  with  (3-38)  and  the  definition  of  ^ 
are  then  substituted  into  the  first-order  continuity  equation  (3-7) •  After 
a  lengthy  algebraic  manipulation,  the  following  equation  is  obtained. 


F  ^ J+<K«)  [(5-Km)OC«  T:/ 


(3-4l) 


^l8  second>order,  linear  differential  equation  for  Uf  vas  not  given  by 
Napolltano  since  his  analysis  was  confined  to  region  B  within  the  frozen 
expansion  fan.  In  this  region,  where  TTp.  ,  the  coefficient  of 

is  zero  and  Bq.  (3*^l)  reduces  to  a  first-order  equation  for 
TJi  (vhich  Indeed  it  must,  since  only  one  boundary  condition,  17,  s  q  , 
is  available).  'By  taking  advantage  of  the  flrozen  flow  relations,  Eq.  (3-4l) 
may  be  further  simplified  to 


vdiere 


(3-i^3) 


Subject  to  the  boundary  condition  -  O  ,  the  solution  to  Bq,  (3-^l) 

may  therefore  be  written  as  the  following  integral  expression  of  known 
functions : 


u  =-3^f£Lv  pf 


■O-^-eSr:  (3JA) 


This  equation  is  equlvedent  to  the  solution  given  by  Napolitano.  In  con- 
tzast,  the  solution  for  Jjf  in  region  C  is  much  more  difficult,  since  the 
second-order  term  does  not  drop  out  Cl  ^  O  for  *6^  s-e-  i-v). 

In  this  case,  we  have  a  two-point  boundary  value  problem  involving  the 
initial  value  of  U,  given  by  Bq.  (3-Wf)  the  condition 
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(3-45) 


=  f, 


since  V,  •  O  at  the  wall.  With  cTp. ,  “T^  and  ^  constant  for 
(3-4i)  may  he  reduced  to 


(3-46) 


where  the  rl£^t  hand  side  is  a  ccmstant.  Ibfortuaately,  a  transfomatlon 
of  Bq.  (3-46)  to  either  normal  or  self-adjoint  fonn  does  not  yield  a  dif¬ 
ferential  eqiiation  that  can  he  solved  in  terms  of  quadratures.  !Iherefore, 
to  obtain  the  flow  field  downstream  of  the  frozen  fan  region,  (3-46)  must 
he  solved  nmerically. 


VflLth  of,  ,  and  ti,  known,  the  remaining  first-order  variables  Vf  , 
f’  ,  oi*  and  are  determined  algebraically  frcn  the  relations 


v;  =  x 


(3-47) 


^  r  -  Op  v.Mi  (  -hVp  V,-!-  ^  %) 


(3-48) 


/ 


■^+Up  W,  ) 


(3-49) 


0^1 


(3-50) 


all  of  idilch  properly  vanish  at  the  initial  expansion  fan  characteristic. 


24 


3-^  Approximate  Solution  for  Hypersonic  Plov 


Before  dlscxisslng  nvnerlcal  examples  of  the  foregoing  exact  solutions.  It 
Is  of  Interest  to  point  out  a  simplified  approximate  solution  for  I 

vhlch  has  been  found  to  bring  out  many  features  of  the  first-order  nonequl- 
llbrlun  flow  field  reasonably  well. 


Since  the  first-order  atom  concentration  may  be  obtained  analyti¬ 

cally  In  closed  form,  It  would  be  of  Interest  to  obtain  approximate  closed 
form  analytical  solutions  for  the  remaining  first-order  thermodynamic 
variables  In  terms  of  .  Qiese  solutions  may  be  obtained  for  M«  >>/ 
by  simplifying  the  analysis  as  follows.  We  may  observe  that  the  present 
problem  Is  analogoxis  In  many  respects  to  the  nonequillbrlu&  expansion  of 
a  dissociated  gas  In  a  hypersonic  nozzle.  Now  in  the  latter  case.  It  has 
been  established  that  large  departures  from  equlllbriun  have  a  relatively 
small  effect  on  both  the  density  and  flow  velocity  In  ccmpeurison  to  the 
corresponding  effects  on  the  local  composition,  temperature  and  pressvure 
as  long  as  the  fraction  of  the  total  energy  involved  In  dissociation  Is 
not  large  (5  25  percent)  (  Ref.  15).  therefore.  In  view  of  the  aforemen¬ 
tioned  analogy,  It  may  be  argvted  that  the  nonequilibrium  effects  on  density 
and  the  resviltant  velocity  In  the  present  problem  may  be  neglected  for  the 
purpose  of  estimating  P,  and  'tT,  when  the  Incoming  stream  Is  hypersonic 
and  Its  dissociation  energy  Is  not  too  large  a  fraction  of  the  total  enereor* 
Then  on  the  basis  of  the  assmptions 

<Ti  ^  O 


we  may  dispense  with  two  of  the  first  order  flow  equations  [conveniently, 
the  differential  equations  (3-7)  and  (3-37)]  and  obtain  the  following 
purely  algebraic  solutions  for  and  ^  : 


(3-52) 


(3-53) 


-  ;2je-  7  oc. 


P 


(3-54) 


These  closed-fonn  approximations  are  questionable  in  region  C;  for  example, 
(3-54)  does  not  correctly  predict  that  is  a  constant  in  this  region 
and,  in  fact,  yields  at  the  vail  (VV*  0  )•  Nevertheless,  as  will 

be  shown  below,  these  formulae  do  bring  out  the  essential  feattires  of 
hypersonic  nonequilibrim  flow  within  region  B  reasonably  well. 

3.5  Nimerical  Results 

In  this  section,  some  nixnerlcal  results  for  the  first-order  noneqxiilibrlvn 
flow  field  will  be  presented.  We  shall  attempt  to  bring  out  the  main 
ihysical  features  of  the  noneqtiillbrlm  effects,  evaluate  the  Influence  of 
the  various  themochemical  parameters,  and  assess  the  accuracy  of  the 
approximate  solution  for  hypersonic  flow  discussed  in  the  foregoing  section . 

In  figs.  2a  throuc^  2f,  angular  distributions  of  the  frozen  and  first- 
order  variables  in  region  B  are  presented  for  typical  values  of  the  physi¬ 
cal  parameters  and  three  different  Mach  nunbers.  Ihese  curves  pertain  to 
all  possible  turning  angles  v”.  Ohe  corresponding  values  of  for  each 
■V"and  [from  Bq.  (3-32)]  are  shown  in  the  Insert  of  Fig.  2a.  It  is 
observed  that  the  sign  of  the  atem  concentration,  temperature,  pressure 
and  vortlclty  perturbations  is  cemmensurate  with  a  recomblnatlon-dcminated 
behavior  throxi#out  region  B  {fi(,  ko  t  %  >  Fi  >  O  >  ^>0  )•  density 
perturbation  is  initially  cempressive  (07  >0  )  but  subsequently  becomes 
negative  throu^out  much  of  the  frozen  fan  region  \siless  is  quite 
small.  In  the  examples  shown,  it  is  clear  that  the  first-order  flow  field 
is  highly  nonlinear  and  rotational  exceprt  in  the  case  of  very  small  total 
turning  angles  (the  range  of  linearity  decreasing  rapidly  with  increasing 


Generally  speaking,  the  perturbations  Increase  very  sharply  at  the 
beginning  of  the  expansion  and  vlU  In  many  cases  attain  a  peak  value  In 
region  B  before  the  last  frozen  fan  characteristic  Is  reached.  Ihe  Initial 
rate  of  growth  Increases  with  M*;  however,  the  maximvBi  values  are  by  com¬ 
parison  rather  Insensitive  to  .  Both  oj  and  Pi  Increase  more 
rapidly  than  oC^  and  TT,  ,  while  ^  rises  gradually  and  lags  considerably 
behind  the  other  perturbations.  For  ^  5>  the  radial  velocity  pertur¬ 
bation  Is  extremely  small  throvi^out  region  B. 

Hie  Influence  of  the  three  parameters  ,  Kp  and  UJ  on  the  first-order 
distribution  functions  Is  shown  In  Figs.  3,  ^  and  5>  respectively.  In 
Fig.  3,  the  noneqvLLllbrlvm  effects  are  seen  to  be  quite  sensitive  to  ; 
a  20  percent  reduction  In  this  parameter  produces  a  ten-fold  Increase  In 
the  magnltide  of  the  perturbations.  Since  the  flow  Is  recombination- 
dominated,  the  Influence  of  Is  prljnarlly  the  result  of  a  change  In 
the  free  stream  dissociation  level  [see  Bq.  (3-3^) ]•  An  Increase  In  , 
as  shown  In  Fig.  4,  also  increases  the  magnitude  of  the  perturbations  for 
the  same  reason.  In  Fig.  5,  an  increase  in  the  recombination  rate  tempera- 
t\xre  exponent  OJ  from  -1.5  to  zero  Is  seen  to  produce  roufjily  a  25  percent 
reduction  in  oc,  and  T,  (with  far  less  of  a  reduction  in  the  remaining 
perturbations).  Diis  effect  follows  from  the  fact  that  the  local  recom¬ 
bination  rate  has  been  reduced.  It  is  clear  that  any  value  of  u)  between 
-2  and  -1  will  yield  essentially  the  same  perturbation  function  solutions; 
hence,  the  major  effect  of  changes  in  the  value  of  U>  lies  In  the  small 
parameter  P  [Bq.  (2-22)]  and  is  proportional  to  1^'^. 

Shown  In  Fig.  6  are  the  distributions  of  TT,  ,  P|  and  ^  based  on  the 
exact  solution  for  oC,  (^)  and  the  apprcKimations  af  o’,  91  Wp  ^  O 

[Bqs.  (3-52)  throu^  (3-5^)].  Ihese  simple  closed-form  approoclmations 
provide  a  satisfactory  qualitative  description  of  the  flow  in  these 
examples,  and  are  in  good  quantitative  agreement  with  the  exact  solxitions 
in  the  vicinity  of  ^  for  all  three  Mach  manbers  shown.  However, 

the  Initial  rise  and  possible  maxima  In  the  perturbations  are  underesti¬ 
mated,  so  that  (3-52)  throu^  (3-5^)  will  be  inaccurate  throujjiout  region 


B  for  expansions  with  very  small  ir-.#  This  result  Is  to  he  expected, 
since  It  Is  seen  from  Pig.  2a  that  a;  s  O  Is  a  poor  approximation  In 
the  Initial  portion  of  region  B. 


PSxact  nimerlcal  solutions  for  the  complete  first-order  nonequillhrlun  flow 
field  In  region  C  downstream  of  the  last  frozen  expansion  fan  characteris¬ 
tic  were  not  carried  out.'***  However,  an  exact,  closed  form  solution  for  OCf 
has  been  given  [Bq.  (3-36)]  and  the  corresponding  value  of  ^  Is  known 
from  the  foregoing  solutions  In  region  B  [Eq.  (3-^)]. 


As  an  Illustrative  example  of  the  complete  first-order  nonequlUbrlvni  flow 
field  In  region  B,  Pigs.  7a  through  7d  present  a  composite  solution  for  a 
typical  case  at  *10  expanded  throu^  the  maximum  turning  angle 
19*2  degrees  with  the  following  free  stream  conditions: 


/  O 


2  0^000 


r  =  3  45(?  *f< 


“po  —  I  A 7TM  —  2 1 1* 7  pyt. 


16  .6 

a  _l9.S>flO  Cm  p  -  £  Ft-ub 

-.-u, 


~<fSoo  •K)*^  ’W.L.-Sac 


r'^O.MZ  ~  2-22) 


nie  general  effect  of  the  nonequillbrlutt  perturbation  on  the  frozen  thermo- 
chemical  field  can,  of  course,  be  anticipated  from  the  foregoing  analytical 
and  nunerlcal  results  (l.e.,  the  marked  change  in  tr ,  oC  and  P  and  the 


« 

In  this  case  a  linearized  theory  such  as  Bef .  10  can  be  used. 

•* 

These  solutions  are  currently  being  programed;  results  will  be  reported 
in  a  future  report. 


less  significant  change  in  <r  ) .  However,  the  graphical  presentation  (Figs . 
7a  through  7d)  of  the  canposlte  flow  field  clearly  shows  to  what  degree 
the  nonequilibrlTia  effects  may  alter  the  basic  field.  In  the  initial  por¬ 
tion  of  the  fan,  the  nonequilibrlvin  effects  indicate  a  flow  field  recom¬ 
pression  f or  r  >  I  inch;  this  recompression  is  unrealistic  and  will  not 
appear  in  the  physical  flow  field.  It  follows  that  in  the  initial  portion 
of  the  fan  the  radius  of  validity  of  the  first-order  solution  is  of  order 
1  or  less  and,  in  fact,  at  larger  radii  the  indicated  composite  solution 
has  been  forced  to  relax  throxigh  ftnd  away  from  the  equilibrium  solution, 
nils  excessive  relaxation  results  in  a  completely  unrealistic  solution. 

As  we  proceed  into  the  fan,  however,  the  radius  of  validity  of  the  solu¬ 
tion  is  extended  outward  since  the  basic  (frozen)  flow  field  departs  more 
and  more  from  equlllbrlm.  In  the  latter  part  of  the  fan  (approaching  the 
wall),  the  nonequllibrltmi  effects  become  negligible  since  the  flow  field 
becomes  frozen  to  all  orders  as  it  approaches  its  maximum  expanded  state. 
Qie  behavior  in  the  initial  portion  of  the  fan  clearly  emphasizes  the 
important  point  that  the  application  of  any  small  perturbation  technique 
stich  as  is  used  in  this  analysis  is  not  simply  a  natter  of  multiplication 
and  addition.  It  is,  rather,  a  careful  process,  foremost  in  which  is  the 
detezmlnatlon  (or  estimation)  of  regions  of  reasonable  validity.  Ihese 
regions  generally  can  be  described  only  within  the  context  of  a  specific 
problem  and  are  dependent  not  only  on  the  Initial  physical  problem  posed 
but  also  on  the  consistency  or  accuracy  of  the  solution  required. 

U.  CONCLUDIIK}  REMARKS 


The  problem  of  the  centered  expansion  of  steady,  supersonic,  chemically 
reacting  diatomic  gas  flow  around  a  sharp  comer  has  been  analyzed  by 
treating  small  local  departures  frm  the  frozen  flow  near  the  comer. 

An  exact  solution  utilizing  a  Idg^hiU-type  gas  model  was  presented  in 
similitude  fom  together  with  some  nunerical  results  for  the  first-order 
nonequlUbrlutt  flow  field.  Ihls  solution  and  the  accompanying  nunerical 
results  brought  out  the  following  features.  Ihe  first-order  species  equa¬ 
tion  was  placed  in  closed  form  and  was  shewn  to  be  independent  of  all  b^tt 
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the  frozen  solution  In  the  expansion  fan  and  at  the  downstream  vail,  acquir¬ 
ing  a  dependence  on  the  previous  first-order  species  solution  In  the  region 
between  the  final  fl*ozen  Characteristic  and  the  wall.  Qie  first-order  var¬ 
iations  of  the  thenao-chemlcal  parameters  show  the  flow  field  to  be  heavily 
reccmblnation-damlnated  but  not  to  the  degree  where  the  neglect  of  the 
dissociation  effects  will  give  an  accurate  assessment  of  the  flow  field. 
The  first-order  flow  field  is  shown  to  be  highly  nonlinear  and  rotational 
except  in  the  case  of  very  small  flow  turning  angles.  Qie  small  range  of 
linearity  decreases  rapidly  with  increasing  free  stream  Mach  number.  In 
general,  the  noneq.ulllbrlum  perturbations  Increase  very  sharply  In  the 
Initial  portion  of  the  expansion  and  will  In  many  cases  obtain  a  peak 
value  before  the  final  frozen  characteristic  Is  reached.  Die  initial  rate 
of  growth  Increases  with  increasing  free  stream  Mach  nmber,  but  the  peak 
values  are  shown  to  be  rather  Insensitive  to  the  free  stream  Mach  nxmber. 

Since  It  has  been  observed  that  for  reacting  expending  flow  in  a  hypersonic 
nozzle  large  departures  from  equilibrium  cause  relatively  small  changes  In 
density  and  flow  velocity,  an  analogy  was  made  with  the  present  problem. 

An  approximate  solution  was  developed  based  on  this  analogy  which  allowed 
a  completely  closed  form  detexmlnation  of  the  thermochemical  field.  These 
simple  closed  foxm  approximations  provide  a  good  qualitative  description 
of  the  flow  field.  Vhen  compared  to  the  exact  first-order  solution.  It  Is 
seen  that  the  Initial  rise  and  the  maximum  values  of  the  perturbations  are 
underestimated.  However,  the  exact  and  approximate  solutions  are  In  good 
quantitative  agreement  In  the  middle  and  latter  portions  of  the  fan. 


DISTRIBUTION  OF  DYNAMIC  AND  THERMODYNAMIC  VARIABLES 
WITHIN  THE  FROZEN  FAN 


TURNING  ANGLE 
V  (DEG) 


DISTRIBUTION  OF  THERMODYNAMIC  VARIABLES 
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DISTRIBUTION  OF  THERMODYNAMIC  VARIABLES 
WITHIN  THE  FROZEN  FAN 


-  1  0 


Kp  .  2.221  10 


DISTRIBUTION  OF  THERMODYNAMIC  VARIABLES 
WITHIN  THE  FROZEN  FAN 
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RGURE  4A 


DISTRIBUTION  OF  THERMODYNAMIC  VARIABLES 
WITHIN  THE  FROZEN  FAN 


-  1  0  Tp  .  17.1 

<u--1.5  A  -  0  K|^  -  3.5 


! 


I 

I  I 


nOURE  4B 


37 


DISTRIBUTION  OF  THERMODYNAMIC  VARIABLES 
WITHIN  THE  FROZEN  FAN 


M.-10  y«  =  1.59  r, 

-  .667  A  -  0  K 
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COMPARISON  OF  EXACT  AND  APPROXIMATE  FIRST  ORDER 
SOLUTIONS  WITHIN  THE  FROZEN  FAN 
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COMPARISON  OF  EXACT  AND  APPROXIMATE  FIRST  ORDER 
SOLUTIONS  WITHIN  THE  FROZEN  FAN 
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FIRST  ORDER  NONEQUILIBRIUM  FIELD 
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ENTROPY  AND  gREE  ENEBOY  OP  A  DISSOCIATED  DIATCKEC  OAS 


The  specific  entropy  of  a  theznally  euid  calorically  perfect  gas  species  Is 

Sj,  S  ^  “pi  ,  (A-l) 

vhere  7e  Is  an  arbitrary  reference  temperature,  Is  the  entropy  at 

To  and  one  atmosphere  pressure,  and  Is  the  partial  pressure  In  atmos¬ 
pheres.  The  corresponding  specific  entropy  for  a  dissociated  diatomic  gas 
mixture  is  thus 

s = (is)  ]  (^) 


vhere  ve  have  taken  "Tm  for  the  expansion  problem.  The  free  energy 
difference  is  obtained  frcm  Eqe.  (2-8),  (A-2)  and  the  definition 


UBi«g 

Now  for  thezmodynamlc  equlUbrivn,  both  the  reaction  rate  function  Gr  in 
Ell.  (2-13)  and  oust  vanish.  Hence  a  comparison  of  B^s.  (2-13) 

and  (A-3)  at  equillbrlUB  diows  that  the  constants  A  and  B  In  Eg.  (2-13) 
may.  In  fact,  be  Identified  with  the  vibrational  specific  heat  parameter 
Kftf  and  the  standard  component  entropies  as  follows: 

As5'-K’„  B  =  ^  -  A  (AJ») 
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Therefore,  with  the  aid  of  Eqa.  (A-4)  and  (2-13),  Bl-  (A-3)  nay  be  re¬ 
written  in  terma  of  G-  : 


r 


G/r 

SirArS')  * 


(A-5) 
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EVALUATION  OF  FIRST-ORDER  DISSOCIATION  IWEQRAL 

Ihe  dissociation  Integral  [Eq.  (3-3519)]  may  be  rewritten  in  a  more  useful 
Ibm  by  eaqploylng  the  trlgononetrle  stibstltutlons  Sgc'y 
Oxen  letting  Wpr-f-2ty-2A  »  Cos*)**  *04  X  * Taw ^ 

obtain  Eq.  (3-35b)  as  follows: 

j  (1*-^ 

It  is  Imnediately  seen  that  When  A/p  ■  2.  ,  the  lotegraxl  (B-l)  reduces  to 
the  ezTor  function;  thus 

••Ap 

x,«,  w  Tan  %)J  (BS) 

IfNp  is  a  positive,  even  Integer  >2,  Bg[.  (B-l)  may  also  be  expressed  in 
closed  form  by  taking  advantage  of  Eq,.  (B-l)  with  NpsZ  and  the  follow¬ 
ing  deilvatlves  thereof: 

-2X/ 

(B-3) 
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e  oiy 


(B-l) 
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!Qius,  for' exanq^le,  ve  have 


VJJTam'V  Tow'll 

vr.r*«l(t 


% 


(B-U) 


-A|> 


^  "  /•  - 

•^pUny^ 


In  this  vay,  one  may  proceed  to  evaluate  Bi.  (B-l)  for  Np  •  6,  6,  ..., 
etc.,  using  progressively  higher  derivatives  frcn  E^.  (B-3)* 
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